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1. INTRODUCTION 

The power shortage relative to the cumulative demand for load is considered to be a major problem 
in many countries where it has become impossible to generate sufficient energy using traditional means. 
These difficult situations have led researchers to focus on finding alternative ways to generate energy [1, 2]. 
Power electronic converters such as DC-DC converters (sometimes called switching regulators) are primarily 
implemented to improve energy conversion efficiency when extracting electric power [3]. DC-DC converters 
are circuits which typically supply a constant output and convert DC voltage to a different voltage level. 
These are used for obtaining stabilized or changing DC voltage(s) by increasing, decreasing, or multiplexing 
from any DC source [4, 5]. 

There are two types of DC-DC converters, isolated and non-isolated. Flyback and forward 
converters are variants of isolated DC-DC converters [6]. These use a high-frequency transformer to place an 
electrical barrier and isolate both the input and output of these converters. A significant advantage of isolated 
converters is that they protect sensitive loads [7]. On the other hand, non-isolated DC-DC converters have no 
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such electrical barrier. Some examples of non-isolated converters are buck-boost converter, buck, boost, 
CUK, ZETA, and SEPIC converters [8, 9]. Non-isolated DC-DC converters are cost-effective and easy to 
design when compared to their isolated counterparts. 

The converters have diverse applications in electric traction, electric vehicles, distributed-DC 
systems like space applications, ships, and airplanes [10, 11]. Solar photovoltaic and specialised electrical 
machine drives are other areas where the converters are useful [12, 13]. Different systems have specific 
requirements and the DC-DC converter used for the applications should be carefully chosen to have a high- 
efficiency system with excellent power quality [14]. This paper focuses on the single-ended primary- 
inductance DC-DC Converter (SEPIC). The SEPIC converter is able to reduce or raise the electrical potential 
(voltage) at the output and can be considered as a buck/boost converter [12]. Additionally, the SEPIC DC-DC 
converter is also capable of producing a regulated positive output voltage for any input voltage. This is in 
contrast to other converters like CUK or ZETA that provide a negative regulated output voltage [15—17]. The 
structure of the paper is specified as follows: introduction, principle of operation of the SEPIC converter, 
state-space modelling for the SEPIC converter, parameters and design components required for the converter, 
MATLAB simulations and modelling the SEPIC converter, genetic algorithm to optimise the PID parameters 
to control the SEPIC converter, results from the simulation, and discussion. 


2. PRINCIPLE OF OPERATION OF SEPIC CONVERTER 

The operating principle for a SEPIC DC-DC converter used in continuous conduction mode (CCM) 
is described in Figure 1. The diagram shows an input voltage source (Vs), a couple of inductors (Li and L2), a 
coupling capacitor (C1) that is connected to the inductors to ensure that a current path exists for the DC flow 
and to also provide insulation [14], power diode (D;), output filter capacitor (C2), and a switching device. 
Usually, the switching device is an insulated gate bipolar transistor (IGBT) labelled (S1). The load resistance 
of the circuit (R) is also shown. 

The SEPIC converter has two operating modes: one is the CCM, the other is the discontinuous 
conduction mode (DCM). In this paper, the first mode is considered to have two states, which can be 
described as follows. 


2.1. SEPIC converter operating modes 
2.1.1.State I (0 < t < D) 

When S1 is switched ON, the current IL; increases. Current flowing through the second inductor l2 
also increases but in the negative direction. The energy stored in the first inductor L; increases. Therefore, 
capacitor Cı supplies energy to increase the current in IL2. SEPIC converter operation mode with switch S; 
ON is shown in Figure 2. 
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Figure 2. SEPIC converter operation mode when S1 is ON 
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The corresponding equations to state I of CCM using KVL and KCL are shown in (1), (2), (3), 








and (4). 
diL1 
Ci —— + lii = 0 (2) 
L —* Va =0 (3) 
dVc2 Vc2 
C2 Ge tR 5 
2.1.2.Stae I (0 < t < D) 


When switch S; is turned off, the capacitor input current (Ici) equals the current through the inductor 
(iri). Current l2 remains in the negative direction and does not reverse its direction. Diode D, actively 
conducts. SEPIC converter operation mode with switch Sı OFF is depicted in Figure 3. 








Figure 3. SEPIC converter operation mode when S4 is OFF 


The corresponding equations to state II of CCM using KVL and KCL are shown in (5), (6), (7), 





and (8). 
Vs = L; 2 Ly Voit Vo2 (5) 
| pee Tag J a sÜ (6) 
a 
Cy a — i =0 (7) 
C2 ez — 1p, — 12 + 2n 0 (8) 


2.2. State space modelling of SEPIC DC-DC converter 


The general state space equations for the SEPIC operating under state I of CCM can be given as in 
equations (9) and (10). 


x(t)= A; x (t) + By u(t) (9) 
y(t) = Cı x(t) (10) 


Then, the state matrix for the SEPIC operating under state I is given in (11) and (12). 





dita 0 0 0 0 
dt 1 š 1 
dir2 0 0 > 0 3 A 
dt —1 L2 0 
aes 0 = 9 O lyn tt i Vs (11) 
dt 1 
dv} | 9 0 O za| œ Lo 
dt 
Und 
Ll 
yit)=[0 0 0 1] a (12) 
Vc2 
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Similarly, the general state equations for the SEPIC converter operating under state II of CCM are 
given in (13) and (14). 


X(t)= A2 x (t) + B2 u(t) (13) 
y(t) = C2 x(t) (14) 


Then, the state matrix for the SEPIC converter operating under state II can be given as shown 
in (15) and (16). 








diy [o 0 = = 
dt z 1 
diz| |0 0 0 = ~ 
dt = L2 0 
dvaj] 0 0 0 jival lo MS (22) 
dt 
avc2 i i 0 Z1 Vez 0 
dt C2 C2 RC2 
lra 
l 
y(t)=[0 0 0 1] a (16) 
Vc2 


The average form of the state matrix over a switching cycle for the SEPIC converter is as given 
in (17) and (18). 








d F 1 
diz} | oq o 2 = |ia] 
dt =) L2 L2 li> 0 
deel ee =R 0 0 Ilva Vs (17) 
i C1 C1 : 0 
aves) || 222 1-D go = {beet Lo 
dt C2 2 RC2 
lr4 
Ll 
y)=[0 0 0 |; (18) 
C1 
Vo2 


Where A represents the system matrix, B represents the input matrix, C represents the output matrix, 
and D represents the duty cycle of the SEPIC converter within the range of (0-1). 


2.3. Component design of SEPIC DC-DC converter 

It is necessary that the selection of the main components in the converter design improve the output 
power efficiency without increasing cost, especially in high-volume power electronic applications. Based on 
the above, this section of paper focuses on choosing components and deriving their mathematical expressions 
when designing a SEPIC DC-DC converter to meet the requirements tabulated in 1. 


2.3.1.Component design of SEPIC DC-DC converter 

The principle of volt-second balance takes into account the voltage drops at diode (Dj), Le. 
according states to the volt-second balance for a converter operating in a steady state, the average voltage 
during one switching cycle of the inductor must be zero [16]. Firstly, the voltage can be seen by Lı during 
interval (0 < t < Dr) equals to (VLi=Vs), while, during the interval (Dr < t < T), it equals (VLi=Vs-Vci-Vp- 
Vo). Taking the average over the entire interval and making it equal to zero, equation (19) can be obtained as: 


1//T 1[ (DT T 
S Viadt) =- fo Vs dt + fo (Vs — Vex — Vp — Vo)dt| = 0 (19) 
Finding the integral and rearranging the above equations yields equation (20) 
V 
Veci + Vp + Vo = — (20) 
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Similarly, by taking the volt- second balance for L2. Hence, the voltage can be seen by Lz during 
interval (0 < t < Dr) is equal to (ViL2=Vc1), while during the interval (Dr < t < T), the voltage is equal to 
Vi2=-VotVp which gives (21). 


1 T 1 DT T 
A(T Viadt) = ESE Ver dt + [* (Vo — Vp) de] = 0 (21) 
Ver + Vo + Vp = 2 (22) 


By combining (20) and (22), the transfer ratio or voltage gain (G) of the SEPIC converter in terms 
of (D) can be obtained as shown in (23). 


G= OME Ds pe (23) 
Vs 1-—D 


The value of (G) can be greater than or less than one, depending on the value of the duty cycle (D). 


2.3.2.Calculation of minimum and maximum duty cycle 
Depending on the voltage gain in equation (23), the minimum duty cycle (D,,;,) occurs when the 
input voltage is at a maximum value (Vs; max) and it can be expressed as shown in (24). 


VotVp 
VotVptVs; max 


(24) 


D pin = 


However, the maximum duty cycle (Dmax) may occur when the input voltage is at a minimum value 
(Vs min), aS given in (25). 


VotVp 
VotVptVs5 min 


(25) 


Dmax ~~ 


2.3.3.Inductor Li and L2 selection 

In switch mode power supply (SMPS), the main function of the inductors is to store energy in their 
magnetic field and attempt to maintain a constant current, or equivalently, to limit the rate of change in 
current flow to the output. Selecting an inductor is necessary to the overall design of a converter. Normally, 
the inductance value of a converter should be set to limit the peak-to-peak ripple current flowing to the 
output. The inductance value that satisfied the requirements in the design of the SEPIC converter was 
estimated by rearranging (1), presented in the previous section, and substituting (di, = Al,, dt = DT, 


T = a to yield (26). 


— 7. — “smin*Pmax 
ly a (26) 


2.3.4. Coupling capacitor Cı selection 

Capacitors have several functions in SMPS design, such as energy storage, filtering, compensation, 
etc. Typically, SMPS stages, the capacitance stores energy as an electric field due to the voltage applied and 
attempts to maintain a constant input and output voltage. To calculate the optimal value of (C1) in the SEPIC 
converter circuit shown in Figure | with an acceptable level of voltage ripple to satisfy the requirements and 
stability of the converter, the capacitance formula shown in (27) was used. 


_ 40 
c= 27) 


Equation (27) can then be rewritten, as the change in voltage is proportional to the change in charge 
over its capacitance, which yields (28). 


_ AQ 
AV = Ta (28) 
However, 
AQ = f i(t)dt (29) 
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Thus, (29) is substituted into (28) to determine the value of input capacitor (C,) as shown in (30). 
Hence, it should be noted that (C,) is charged by (l2 = Iomax) during the time interval (0 <t < DT). 


f . 27 dt 
C = E (30) 
Solving the integration shown in (30) gives (31). 
me !omax*Dmax (31) 


AVc1* fs 


2.3.5. Output capacitor C2 selection 
Likewise, as in the coupling capacitor (C,) calculation, the value of output capacitor (Cz) as 
obtained in (32) should supply the output current (I, max) which is loaded during the ON state. 


ĉ; = Iomax*Pmax (32) 


AVc2*fs 
Where Dmin, Dmax are the maximum and minimum duty cycles respectively , V; max Vs min are 
maximum and minimum input voltages respectively, Vg is the output voltage in volts, and Al, is the same 
peak-to-peak ripple current in amperes for both inductors L, and L, when the minimum input voltage (V; min) 
and switching frequency (f>) are applied, it may be expressed as approximately (AI,=20% *lin max). 4Vc1; 


AV; are the peak-to-peak voltage ripples in volts at capacitors C} and Cj, respectively. The acceptable 
percentage value of ripple voltage used in this paper is (AVc1= 1% *V> min and AVc2=1% *Vo). 


2.4. Matlab simulation of SEPIC DC-DC converter 

The open loop SEPIC DC-DC converter was implemented using the MATLAB® simulation 
platform R2018b Simulink [18]. For the purpose of simulation, the components of the SEPIC converter and 
their values tabulated in the Table 1 depend on the equations derived in pervious sections. 


Tablel. Component values of the designed SEPIC Converter 


No Parameter Symbol Value 
1 Minimum input voltage Vs min 100 V 
2 Maximum input voltage Vs max 120 V 
3 Output voltage Vo 311 V 
4 Maximum output power Po max 2 Kw 
5 Maximum load resistance Rypmax ~50Q 
6 Maximum load current lömaz 6.43 A 
7 Switching frequency fs 25 kHz 
8 Diode drop voltage Va 0.7 V 
9 Minimum duty cycle D min 0.722 

10 Maximum duty cycle Dax 0.757 

11 Current ripple at L4 , Ly Al, 4.01 A 

12 Voltage ripple at C4 AVc1ı 1V 

13 Voltage ripple at C, AVc2 3.11 V 

14 Coupling capacitor Ci 195 uF 

15 Output capacitor C3 63 uF 

16 Inductor 1 Li 755 uH 

17 Inductor 2 Lo 755 uH 


The output to the input voltage transfer function, G(s) can be derived from the average state (17) and 


(18) by using the formula to convert the state-space matrix to a transfer function, shown in (33), and 
substituting the values of the design parameters and components to obtain (34). 


G(s) = C (SI-A)! B (33) 
Vo (S) 5.135*108 S? + 1.824*1076 S+ 2.642*1015 

a 4 3 2 9 12 (34) 
Vin(s) S4 +330.2* S3 + 6.792*10© S? + 1.419*109 S + 8.478*10 
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2.5. Genetic algorithm (GA) 

GAs are stochastic global search engines that simulate the creation of natural processes. A GA starts 
with an initial population of chromosomes, each of which represents a possible solution to the problem, and 
the reliability of the chromosomes is evaluated by a fitness function (FF) [19, 20]. The GA consists of three 
main steps: mutation, crossover and selection [21]. GAs are used to determine the representation of the 
chromosomes that are created by three values corresponding to the gains or parameters of the PID controller 
(Kp, Ki, Ka ) that have been modified to attain a desired behavior [22]-[24]. The (Kp, Ki, Ka ) gains are real 
numbers, which are to be measured individually. The equation for the PID controller is shown as (35) [25]. In 
this study, the integral time absolute error (ITAE) is taken as the objective function, given as (36). The 
reciprocal of the objective function is called FF. The FF is the measurement of the chromosome’s quality. 


G.(s) =K + +K4s (35) 
ITAE =f t l e(t)ldt (36) 


Where K, is the proportional gain, K; is the integral time, Kg is the derivative time, e (t) is the error- 


controlling signal = 1— y(t) and y(t) is the tuned control system step response. Figure 4 shows a flowchart for 
GA. The parameters of the GA are listed in Table 2. 


Table 2 Genetic algorithm parameters 












Start 
No Parameter Type / Value 
1 No. of iteration 100 
2 Population size 50 
3 Encoding Binary Create initial population as 100 
4 Selection scheme Stochastic Uniform Chromosome 
5 Crossover probability (Po) 0.8 initial value of [Kp, Ki, Kd] 
6 Crossover type Single Point 
7 Mutation probability (P,,) 0.2 | 
8 Recombination probability (P,) 0.09 Evaluate the fitness for cach 
9 Fittness function ITAE Chromosome 


| 


Perform Selection, Crossover, 
and Mutation Process 





JFitness $ € 


| Yes 


Best Chromosome 
[Kp, Ki, Kd] 


| 


END 


Figure 4. Flowchart for genetic algorithm (GA) 


The block diagram for the GA-optimized or tuned PID controller parameters for the SEPIC DC-DC 
converter can be seen in Figure 5. The reference signal or set point R(s) is input as the desired voltage and is 
compared with the actual value of the output voltage Y(s) for the SEPIC converter. The error signal E(s) is 
given as the input signal to the PID controller after its parameters are tuned using GA technique. The output 
signal generated by the PID controller U(s) is used as the input signal to control the suitable range of duty 
cycle (D) in a suitable range for the SEPIC converter, and the required output voltage is then obtained [26]. 
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Figure 5. Block diagram of PID controller optimized by GA 


3. RESULTS AND DISCUSSION 

MATLAB® R2018b Simulink software was used to simulate the SEPIC DC-DC converter and 
verify its performance. Figure 6 (a) shows the step response of the open-loop transfer function for the SEPIC 
DC-DC converter without the PID controller. It is observed that the response matches the step response for a 
second-order system. However, the maximum overshoot percentage (close to 73%) is higher than that of the 
second-order system. Slow rise and settling times are also observed. Figure 6 (b) depicts the closed-loop step 
response for a GA optimised PID controller. The PID controller parameters (Kp, Ki, Ka ) are set to their 
optimal value using the GA technique. This technique helps minimise the objective function (ITAE) and the 
optimisation process depends on the code specified in the MATLAB m-file. Table 3 specifies the 
optimisation criteria and the performance characteristics for the step response like (Mp, tr, ts). 


Step Response of Opened Loop TF for SEPIC Converter Step Response for Closed Loop with Optimized PID Controller Based GA 


| GA-PID | 





Without PID Control ler | 





Amplitude 
Amplitude 


Time (seconds) Time (seconds } 


(a) (b) 


Figure 6. Simulation results, (a) Step response for SEPIC converter without controller, (b) Step response for 
SEPIC converter with GA-PID controller 


Table 3 Performance Parameters for step response 


Performance Type of controller 
Parameters Without PID controller With GA-PID controller 

K, N/A 69.805 
K; N/A 290.779 
Ka N/A 339.306 

% Mp 73 0 

t, (sec) 8.9246 e-04 1.2610 e-11 

t, (sec) 0.0301 2.2454 e-11 

ITAE N/A 2.7182 e-06 
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4. CONCLUSION 

This paper presents design, modelling and simulation of controlled SEPIC DC-DC converter-based 
GA. The study comprised the mathematical model, state-space modelling, design components, simulation, 
and analysis of the SEPIC DC-DC converter. MATLAB R2018b SIMULINK was used to carry out the 
simulations. The SEPIC DC-DC converter is a time-invariant, linear system of the fourth order. A PID 
controller was utilised to change the output of the SEPIC converter to optimise the values of the performance 
parameters. PID controller parameters Kp, Ki, Ka were manipulated using the GA technique to determine 
their optimal values intelligently. Optimisation was performed with ITAE serving as the fitness function. The 
results of the simulation and the study indicate that the closed-loop step response system performed best in 
conjunction with a PID controller with optimised parameters. Maximum overshoot percentage rise and 
settling times were optimised. Additionally, the GA-PID controller was more accurate and better sensitivity 
against disturbances when compared to using it without the PID controller. 
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